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ABSTRACT 


The coalescence of massive black hole (BH) binaries due to galaxy mergers 
provides a primary source of low-frequency gravitational radiation detectable by 
pulsar timing measurements and by the proposed LISA (Laser Interferometry 
Space Antenna) observatory. We compute the expected gravitational radiation 
signal from sources at all redshifts by combining the predicted merger rate of 
galactic halos with recent measurements of the relation between BH mass, Mbh, 
and the velocity dispersion of its host galaxy, a. Our main hndings are as follows: 
(i) the nHz frequency background is dominated by BH binaries at redshifts z ^2, 
and existing limits from pulsar timing data place tight constraints on the allowed 
normalization and power-law slope of the Mbh-cr relation or on the fraction of BH 
binaries that proceed to coalescence; (ii) more than half of all discrete mHz mas¬ 
sive BH sources detectable by LISA are likely to originate at redshifts z >7; (iii) 
the number of LISA sources per unit redshift per year should drop substantially 
after reionization as long as BH formation is triggered by gas cooling in galaxies. 
Studies of the highest redshift sources among the few hundred detectable events 
per year, will provide unique information about the physics and history of black 
hole growth in galaxies. 
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universe 
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1. Introduction 


A future detection of gravitational waves (GWs) will not only test the theory of Gen¬ 
eral Relativity, but will also provide a precious new tool in astronomy (Hughes et ah 2001, 
and references therein). This tool will be particularly effective for studying compact ob¬ 
jects because the information carried by GWs is very different from that of electromagnetic 
radiation. Gravitational radiation is generated by different physical processes and is often 
emitted from dense regions that are otherwise hidden from traditional astronomical obser¬ 
vations. A very important distinction, particularly for the study of the distant universe is 
that while the usual observable for electromagnetic radiation is an energy flux which drops 
in proportion to distance squared, the direct observable for gravitational radiation is the 
wave amplitude h which drops off only as distance. As a result, the coalescence of many 
massive black-hole (BH) binaries (Mbh > lO^M©) will be detectable by the planned GW 
observatory LISA (Laser Interferometry Space Antenna^) at mHz frequencies out to very 
high redshifts (Haehnelt 1994; Hughes et al. 2001). At the lower /xHz frequencies, pulsar 
timing measurements limit the level of the cosmic GW background, contributed mainly by 
the early phase in the coalescence of massive BH binaries (Kaspi et al. 1994; Thorsett & 
Dewey 1996; Lommen 2002). 

Massive BH binaries are among the primary sources for LISA. Such binaries are produced 
by mergers of galaxies with pre-existing black holes, provided that coalescence due to the 
interaction of the binary with its surrounding galaxy occurs in less than a Hubble time. 
Hughes (2002) has shown that LISA observations of the gravitational radiation during both 
the inspiral and ring-down phases of massive BH binary coalescence can be used to determine 
the source redshift, while the source position on the sky can be inferred within a fraction of 
an orbit of the LISA satellite. Recently, Seto (2002) has noted that consideration of the hnite 
length of the LISA detector leads to signihcant improvements in estimations of the accuracy 
to which the three-dimensional positions for coalescing BH binaries can be measured. 

Black holes are ubiquitous in low redshift galaxies (e.g. Magorrian et al. 1998). Their 
masses, Mbh, as determined from dynamical measurements (e.g. Merritt & Ferrarese 2001; 
Tremaine et al. 2002) or from reverberation mapping of AGN (e.g. Ferrarese et al. 2001; 
Laor 2001), correlate tightly with the velocity dispersion, a*, of their host stellar bulge. 
However, the abundance of massive BHs in galaxies at high redshifts is not well known. In 
fact, Menou, Haiman & Narayanan (2001) have argued that an occupation fraction as low 
as 0.01 at redshift 5 can result through subsequent mergers, in all local massive galaxies 
containing massive BHs. Nevertheless, one might expect massive BHs to be ubiquitous at 
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high redshifts for the following reason. It is commonly thought that present day BHs are the 
dormant remnants of the active galactic nuclei observed at high redshift. Recent observations 
of this population have reached out beyond redshift 6 (Fan et ah 2001). The quasars discov¬ 
ered are very bright, and must be powered by BH masses of > IO^Mq. Recent modeling by 
Wyithe & Loeb (2002) has shown that the density of these high-redshift quasars as well as 
lower redshift populations can be explained through merger activity using the same relation 
between host velocity dispersion and BH mass as observed locally. The latter assumption is 
supported by the recent observation that the relation between velocity dispersion and BH 
mass holds in quasars out to redshift 3 (Shields et ah 2002). This result suggests that if BHs 
are not ubiquitous in galaxies at high redshift, then it must be a coincidence that the active 
quasar period increases at high redshift in inverse proportion to the occupation fraction. 
Note also that the inferred active period for a BH at ~ 6 is ~ lO^yr or ~ 1% of a Hubble 
time (Wyithe & Loeb 2002). Thus the minimum occupation fraction at ~ 6 is ~ 0.01. 

The emission of GWs requires that the binary coalesce in less than a Hubble time. What 
is the mechanism that brings a wide binary system with an orbital separation > 1 kpc into the 
regime where coalescence is dominated by gravitational radiation? Several mechanisms have 
been discussed in the literature, including dynamical friction on the background stars (e.g. 
Begelman, Blandford, & Rees 1980; Rajagopal & Romani 1995; Milosavljevic & Merritt 2001; 
Yu 2002), and a gas dynamical effect analogous to that thought responsible for planetary 
migration (Gould & Rix 2000; Armitage & Natarajan 2002). Yu (2002) has followed the 
expected evolution of massive central BH binaries in a sample of nearby early-type galaxies 
for which high-resolution photometry exists, and found that the evolution timescale increases 
signihcantly when the BH binary becomes hard before entering the GW regime. Spherical 
galaxies often lead to coalescence timescales that are longer than a Hubble time, while 
highly flattened or triaxial galaxies lead to faster coalescence. Regardless of the physical 
mechanism, Haehnelt & Kauffmann (2002) argue that if central binaries do not merge in 
less than a Hubble time, then three-body interactions resulting from subsequent mergers 
will cause ejection and too much scatter in the Mbh — cr* relation. Hereafter, we assign 
an efficiency emrg to the coalescence process. Throughout this paper we show numerical 
results under the assumption that BH binaries formed through the merger of two galaxies 
proceed quickly to coalescence (e^rg = 1), but point out the relevant scaling with emrg where 
appropriate. 

The number of binary sources depends on a convolution over redshift of the merger rate 
of galaxies with the fraction of galaxies that harbor massive BHs. Several attempts have 
been made to compute the event rate (Fukushige, Ebisuzaki & Makino 1992; Haehnelt 1994; 
Rajagopal & Romani 1995; Haehnelt 1998; Menou, Haiman & Narayanan 2002; Jaffe & 
Backer 2002), either based on observational counts of quasars and elliptical galaxies with 



arbitrary extrapolations to higher redshift, or using merger tree algorithms based the Press- 
Schechter theory. Despite the fact that detectors such as LISA are expected to detect events 
at exceedingly high redshift, many calculations have only considered low (z < 5) redshift 
sources. However, recent evidence of bright quasar activity at z ^ 6 (Fan et al. 2001) and 
independent evidence that these quasars reside in massive dark-matter halos (Barkana & 
Loeb 2002), suggests that BHs might be as ubiquitous at the centers of high-redshift galaxies 
as they are today. Furthermore, prior to reionization the ability of gas to cool inside small 
halos is mainly limited by the efficiency of gas cooling rather than by the infall of gas from a 
heated and photoionized intergalactic medium, and so massive BHs could form inside much 
smaller galaxies than are available after reionization (e.g. Barkana & Loeb 2000; Haiman, 
Abel & Rees 2001; Thoul & Weinberg 1996). As a result it is possible that a substantial 
fraction of the events detectable by LISA have originated prior to reionization. Our paper 
explores this possibility. 

Recently, Jaffe & Backer (2002) have computed the spectrum of characteristic strain 
resulting from massive BH mergers in the nHz-juHz regime using a phenomenological ap¬ 
proach for the merger rate and BH mass function up to some maximum redshift. They 
found that the logarithmic spectral slope of the stochastic cosmic background of gravita¬ 
tional radiation that results from inspiraling BHs is -2/3, independent of the details of their 
model for the merger rate and BH mass function. This provides a numerical conhrmation 
of the general result derived by Phinney (2001). They also found that the amplitude of the 
spectrum is close to the limit probed by recent pulsar timing experiments, but is sensitive to 
the model assumed for the BH mass function and merger rate. In this paper we make a full 
semi-analytic calculation of the entire merger rate history starting from the highest redshifts 
when the first galaxies are believed to have formed in the universe (Barkana & Loeb 2001, 
and references therein). We employ the formalisms of Press & Schechter (1974) and Lacey 
& Cole (1993), and compute the resulting redshift-dependent BH mass function (including 
the effect of reionization) using the observed relation between BH mass and halo circular 
velocity. We compare our results with current limits on the cosmic GW background at low 
frequencies, and extend the calculation upwards in frequency to the LISA band. 

Throughout the paper we assume density parameter values of = 0.35 in matter, 
Qf, = 0.052 in baryons, Da = 0.65 in a cosmological constant, and a Hubble constant of IIq = 
65 kms“^Mpc“^ (or equivalently h = 0.65). For calculations of the Press-Schechter (1974) 
mass function (with the modihcation of Sheth & Tormen 1999) we assume a primordial 
power-spectrum with a power-law index n = 1 and the htting formula to the exact transfer 
function of Cold Dark Matter, given by Bardeen et al. (1986). Unless otherwise noted we 
adopt a present-day rms amplitude of as = 0.87 for mass density fluctuations in a sphere of 
radius 8/i“^Mpc. 
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The paper is organized as follows. In § 2 we describe calculation of the evolution of 
the BH mass-function following reionization. In § 4, § 5 and § 6 we use the results of this 
calculation to predict the merger rate and hence the detection rate of GW sources by LISA, 
the distribution of event durations, and the resulting spectrum of the characteristic strain. 
Finally, we present our conclusions in § 7. 


2. Evolution of the Mass-Function of Halos Hosting Massive Black Holes 


During hierarchal galaxy formation, an average property of the galaxy population, such 
as the density of galaxies evolves in redshift due to newly collapsing halos, mergers of halos 
and accretion. In this section we discuss the calculation of the rate of newly collapsing 
halos, and then demonstrate a method to calculate the continuous evolution of the massive 
BH mass function. 

Let be the Press-Schechter mass function (number of halos with mass between 

M and M -|- dM per co-moving Mpc) of dark-matter halos at a redshift ; 2 . We include the 
extension of Sheth & Tormen (1999) which improves the agreement with direct numerical 
simulations (Jenkins et al. 2001). The rate of change (with respect to redshift) of the density 
of dark matter halos with masses between Mi and Mi -|- dMi is therefore Let us also 

dehne lin, the number of mergers per unit time of halos of mass AM with halos of 

mass M (forming new halos of mass Mi = M + AM) at redshift (Lacey & Cole 1993). 

We consider mergers of halos having masses greater than Mmin = IOOMq. Addition 
of mass from halos smaller than M^in is treated as smooth accretion. In what follows this 
contribution to the growth of halos must be included to maintain a self-consistent calculation. 
We therefore hnd the component of the change in the density of halos with masses between 
Ml and Mi -|- dMi that is due to accretion rather than to mergers. This is given by the 
overall change in the density of halos with mass Mi, minus the change due to merger of halos 
to form a new halo of mass Mi, plus the change in density of halos of mass Mi that merge 
with other halos (a change that was compensated for in the net The net change in 

density of halos of mass Mi through accretion is then. 


<Pn^c 
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We are interested in the evolution of the mass function of dark-matter halos that contain 
a massive BH. Prior to reionization, the minimum mass halo inside of which gas can cool 
is limited by the cooling rate, and corresponds to a virial temperature above Tmin lO^K 
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if the gas cools through transitions in atomic hydrogen, or above Tmin 10^'^K if the gas 
cools through transitions in molecular hydrogen (Haiman, Abel & Rees 2000). Following 
reionization, the IGM is heated to ~ 10"^K. The Jeans mass is increased by several orders 
of magnitude, and numerical simulations hnd that gas infall is suppressed in larger halos. 
There is some uncertainty about the precise value of the halo circular velocity below which 
gas infall is suppressed (e.g. Thoul & Weinberg 1996; Kitayama & Ikeuchi 2000; Quinn, 
Katz & Efstathiou 1996, Weinberg, Hernquist & Katz 1997; Navarro & Steinmetz 1997). 
We assume that following reionization, gas can cool inside newly collapsing halos if their 
virial temperatures are higher than Tmin = lO^K. We denote by Fkcc the fraction of halos 
which contain a massive BH after crossing the cooling threshhold through accretion. 


Halos forming via mergers are assumed to contain a massive BH if one of the subunits 
used to make them already contained a BH. In addition, if the merger product has a mass 
larger than M( Tm,n ), then a BH could form out of the gas that cools inside the merger 
product. We assume that this happens in a fraction FArg of cases where neither the initial 
nor the accreted halo contained a BH. The mass function of halos containing a massive BH 
therefore evolves according to the following differential equation; 
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( 2 ) 


where 0 is the Heaviside step function. If Face = FArg = 1, this equation maintains a density 
of BHs in halos larger than M(T m,r, ) that is equal to the Press-Schechter density. In the above 
equation, the hrst term gives the change in the density of halos containing BHs due to newly 
collapsing halos above the threshold virial temperature. The second and third terms give the 
increase due to mergers of the density of halos of mass Mi that contain a BH. The second 
term gives the increase for halos below the threshold temperature due to mergers with a halo 
already containing a BH. This term will be important since the mass corresponding to the 
minimum temperature increases with time, and because reionization results in an increase 
of the minimum mass halo inside which gas can collapse and cool. The third term gives the 
increase in the density of halos with virial temperatures above the threshold temperature 
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that contain BHs. There are two parts to this term. If the halo of mass Mi — AM contains 
a black hole, then so does the merger product. However if the halo of mass Mi — AM does 
not contain a BH then the merger product only contains a BH if one of the following two 
conditions hold: either the accreted halo (of mass AM) already contained a BH, or if it did 
not then a BH forms in a fraction of merger products. Finally, the fourth term allows for 
the loss in the density of halos of mass Mi that contain black holes due to mergers forming 
a larger halo. There are two limiting cases; Face = -^mrg = 1 in which case all galaxies with 
masses above the cooling threshhold contain black-holes, and Face = -^mrg = 0 in which case 
the number of massive BHs drops with time through galaxy mergers. Below we discuss these 
cases in more detail. 

We consider 4 examples for the evolution of These will be referred to throughout 
the paper and are labeled I — IV. In each case the initial values for the mass function of 
halos containing BHs is chosen to be 

• (I) Begin at 2 : = 20 with BHs in all halos having T^ir > lO^K, and assume Fin, = Face = 
Fmrg = T Reionization is assumed to be at 2:reion = 7. 

• (H) Begin at = 20 with BHs in all halos having Tyir > 10"^K, and assume Fjni = 
Face = Fmrg = 1- Reionization is assumed to be at 2:reion = 12. 

• (HI) Begin at 2 : = 20 with BHs in all halos having Tvi^ > 10^'^K, and assume Fjni = 1 
and Face Rmrg 0- 

• (IV) Begin at 2 : = 20 with BHs in all halos having T^i,. > 10"^K, and assume Fini = 
Face = Rmrg = 1- Reiouizatiou is assumed to be at 2:reion = 12, but in difference from 
case H , we use ag = 0.6 rather than the hducial value of ag = 0.87. 

Cases I, H and IV refer to situations where cooling in neutral gas is limited by atomic 
hydrogen, and in which BH formation is ongoing. Case HI refers to situations where cooling 
in neutral gas is limited by molecular hydrogen, and in which the BH population arises from 
“seed BHs” at z = 20. 

We show the four examples of evolution of in Figure 1. Here the density of halos of 
mass M that contain a BH is plotted as a function of redshift (solid lines). The dashed lines 
show the corresponding Press-Schechter mass function, and each pair of curves is labeled 
by the logarithm of the mass in Mq on the left of the plot. The thick dashed line is the 
density of halos with the minimum virial temperature into which gas can collapse, and has 



logio(dn;,h/dM) log,Q(dnj^ydM) 



Fig. 1.— Solid lines: The mass function of halos that contain BHs as a function of redshift for various 
masses (labeled by \ogif^{M/M q)). Dashed lines: The Press-Schechter mass function of halos as a function 
of redshift for various masses. The panels show cases I-IV whose primary parameters are listed at the top 
of each plot. The thick dashed line shows the density corresponding to Tmin as a function of redshift. 
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a discontinuity at Zreion- In cases I and II we see that following reionization, the density of 
halos below the critical virial temperature that contain BHs drops as these halos merge to 
form more massive halos. By redshift 2 ( 2 :reion = 7) and 5 ( 2 :reion = 12), there are no halos 
containing BHs left with temperatures below Tmin- 


2.1. The Black-Hole Mass Function 

The calculation of the BH mass function from the mass-function of halos that contain 
BHs requires specihcation of a relation between the BH and halo mass. Ferrarese (2002) has 
found a tight, power-law correlation in nearby galaxies between the mass of a dark-matter 
halo and the mass of its central BH. Following Haehnelt, Natarajan & Rees (1998) and 
Wyithe & Loeb (2002) we generalise to higher redshifts assuming a BH mass that scales as 
a powerlaw with circular velocity Mbh oc v^- Using the expression for halo mass M as a 
function of Vc (Barkana & Loeb 2001) we hnd 



where Ac = IStt^ -|- 82d — 39(P, d = — 1 and = Hm(l + 2 :)^/[Hm(l + + ^a]- In this 

notation, the relation found by Ferrarese (2002) for BHs in the local universe has eo = 10 “® °® 
and 7 = 5.46. In the remainder of the paper masses with the subscript bh a subscript refer 
to BH masses, while those without refer to halo masses. The merger rates are computed in 
terms of halo masses, but many expressions implicitly contain conversion from BH to halo 
mass or vice versa using equation (4). 

Simple considerations of the regulation of BH growth by feedback from quasar activity 
result in 7 = 5 (Haehnelt, Natarajan & Rees 1998; Silk & Rees 1998). Supporting this 
argument is the recent hnding of Shields et ah ( 2002 ) that the Mbh — <7* relation derived 
from QSOs shows no redshift dependence. The growth of BHs only through hierarchical 
mergers imply 7 = 3. A value of 7 > 3 is therefore required in order to allow BH growth 
through accretion that must accompany quasar activity. With 7 = 5, theoretical luminosity 
functions require eo = 10“®'"^ (Wyithe & Loeb 2002), and these are the values used throughout 
this work unless specihed otherwise. The mass function of BHs is therefore 

dubh ^ dnbh 3 M 
dMbh dM 7 Mbh 

The evolution of is shown in each of cases I-IV in Figure 2. The mass-function 
of halos that contain BHs (^]^) is insensitive to the BH-halo mass relation. However, in 
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z . =7 T . =10000K F=1.0 a„ = 0.87 
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Fig. 2.— The mass functions of massive BHs at various redshifts. The panels show cases TIV whose 
primary parameters are listed at the top of each plot. 
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the examples shown, BHs exist in halos as small as M < lO^M©. Extrapolation of the 
^bh ~ Vc relation down to these small galaxies is not justihed by observations; moreover, it 
yields BH masses comparable to those expected from stars where the formation physics is 
very different. We therefore only show the mass function down to Mbh = IO^Mq. We will 
show that the low frequency GW signal is dominated by BHs in galaxies with M > IO^^Mq; 
however BHs in smaller galaxies are important as seeds for massive BHs in larger galaxies 
at later times. 

The BH mass functions can be compared with those presented in Figure 5 of Kauffmann 
& Haehnelt (2000). The density of the most massive BHs increases monotonically with time, 
while in cases I, H and IV there is little evolution in the density of lower mass BHs before 
reionization. However following reionization, there is rapid evolution at low masses since 
newly formed low-mass halos do not produce BHs, and the density of low mass halos that 
contain BHs drops as they merge to form larger halos. Case HI shows a different evolution; 
since there is no ongoing BH formation the reionization redshift does not leave any signature. 
Case IV with a smaller value of ag = 0.6 yields lower BH densities at all redshifts. 


3. Gravitational Waves from Inspiraling Binary Black-Holes 


As described earlier, we assume that in a fraction emrg of galaxy mergers the resulting 
binary BHs evolve rapidly into the regime where energy loss is dominated by gravitational 
radiation. Dehning this time as fstart, the subsequent evolution of the binary is then described 
by equation (6). It can be shown (e.g. Shapiro & Teukolsky 1983) that the evolution in the 
binaries’ (circular) orbital frequency / with time t obeys the relation 




start 


5 C® (Mbh + AMbh)^'^^ , . _8/3 / .-8/3 f-8/3^ 

256 GV3 MbhAMbh ^ ’ V ^ J ’ 


( 6 ) 


where /start is the orbital frequency at time tstart- Throughout the paper we assume circular 
orbits for which the observed gravitational radiation has a characteristic frequency 


fc 


f 


1 + ;^’ 


(7) 


and characteristic strain, averaged over orientations and polarizations (Thorne 1987) of 


he 


/2y/' G^AMbhAMbh 2/3 

Vlhy' c4i?(;^)(Mbh + AMbh)V3^ ""71 , 


( 8 ) 


where R{z) is the co-moving coordinate distance to redshift ; 2 . The choice of R{z) as the 
distance measure ensures that the energy flux, which is proportional to (hc/c)^, declines in 
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proportion to the inverse square of the luminosity distance. As a binary loses energy, the 
frequency as well as amplitude of the gravitational radiation increase with time. Following 
Hughes (2002) we assume that the inspiral attains its maximum frequency when the two BHs 
are separated by 3 Schwarzschild radii, i.e. /max = c^/(27r6^/^GM)Hz for a Schwarzschild 
BH, and stop our calculation at this point. Hughes & Blandford (2002) find that following 
a major merger, the remnant is rarely rotating rapidly unless the mass ratio is ~ 1. In 
this work we neglect the subsequent higher frequency radiation from the coalescence and 
ringdown phases. 


4. The Merger Rate and the Counts of Gravitational Wave Sources 


The proposed GW observatory LISA will be able to detect many mergers of massive 
BHs out to very high redshifts (Haehnelt 1994; Hughes et al. 2001). It is therefore inter¬ 
esting to calculate the number of detectable mergers expected per unit redshift per unit 
(observed) time on the sky. Previous attempts at this calculation have used empirical local 
rates and extrapolation (e.g. Backer & Jaffe 2002), or have used a merger tree algorithm. 
Menou, Haiman & Narayanan (2001) used a merger tree algorithm to predict the merger 
rate out to z ^ 5, beyond which they found that the Press-Schechter mass function was no 
longer accurately recovered. We have circumvented this problem by integrating the BH mass 
function forward in time, while exactly preserving the Press-Schechter mass function. This 
allows us to explore the merger rate at higher redshifts and identify the possible signature 
of reionization in GW observations. 


The number of detectable BH mergers per observed year per unit redshift across the 
entire sky is 


d^N, 


r-M 


gw 


dtdz 


dM / dAM 0(M, AM, /c, K, z) X S{z, Mbh, AMbh) 


( 9 ) 
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where V is the co-moving volume and G is solid angle on the sky. Note that there must be 
BHs present in both galaxies in order for gravitational radiation to be produced (this being 
the origin of the term in square brackets). The function 0(M, AM, /c, he, z) has a value of 
1 if the merger produces a detectable strain in the LISA band, but is zero otherwise. In 
our calculations throughout this paper we set emrg = 1- Since emrg is our parameterization 
of the efficiency with which dynamical friction brings binary BHs into the GW regime, 
this corresponds to the assumption that all BH binaries coalesce. However, more generally. 


dtdz 
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Fig. 3.— The rate of GW events per unit redshift, per observed year. The panels show cases I-IV whose 
primary parameters are listed at the top of each plot. We assume eo = and 7 = 5. For comparison 

the dashed line shows the event rate where Tmin is assumed not to rise following reionization. 
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At a frequency of fc = 10“^Hz, the expected value of the threshhold sensitivity is 
he ~ 10“^^ for a signal to noise ratio of 5 (needed for conhdent detection of a binary with 
unknown frequency and direction) and one year of observation. We hnd that the typical 
event duration at frequencies between and fc = 10“^'^Hz ranges between a 

few months to a year. We therefore assume detectability if the strain amplitude he > 10“^^ 
at an observed frequency of fc = 10“^Hz. Note that a full calculation of the signal to noise 
ratio of an observation requires comparison of the coherently folded signal power (including 
the evolution in its frequency) to the noise power accumulated over the bandwidth of the 
measurement. By working with the quantity he (Thorne 1987), we approximately account 
for the detection threshold. 


The function S is included in equation (9) because the merger of dark matter halos 
does not imply the immediate merger of the stellar systems within. Colpi, Mayer & Gov- 
ernato (1999) have studied the evolution of satellites in virialized halos. They hnd that the 
orbital decay time is 
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ecay 


= 1.2 ( ^ 
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M +AM 
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where rvir is the initial radius and e ~ 0.5 is the circularity of the of the orbit (Ghigna 
et ah 1998). Note that even though e ~ 0.5 for the initial orbit of the secondary galaxy 
in the dark matter halo, eventually the eccentricity decays rapidly during the GW driven 
inspiral of the binary BH system (e.g. Peters 1964; Hughes & Blandford 2002; Yu 2002). 
The dynamical time may be written from the expressions in Barkana & Loeb (2001) as 
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and amounts to about a tenth of the Hubble time. In order to have the satellite sink in 
less than a Hubble time, the merging galaxies must be of comparable mass. The dynam¬ 
ical evolution of a massive binary BH may be disrupted if there is a subsequent major 
merger before coalescence of the initial binary. This typically leads to the ejection of the 
lightest of the three black-holes (e.g. Haehnelt & Kauffmann 2002; Volonteri, Haardt & 
Madau 2002). To account for this phenomenon in a simple way, we assume that at most 
one coalescence can occur during the decay time fdecay or within a Hubble time, whichever is 
the larger. In mergers where the satellite does not sink within a Hubble time, we assume no 
coalescence. Furthermore, following Kauffmann & Haehnelt (2000) we suppress gas accre¬ 
tion and therefore BH mergers within galaxies with Vc > 600km/sec, though this restriction 
does not signihcantly affect our results. We therefore include the function 5'(2:, Mbh, AMbh) 


which is set to 


max 


(ami 


<fp 

decay dAMdt 


M 


-1 


if the satellite galaxy’s orbital decay time 


fdecay is Smaller than the Hubble time, ( i7oA/Gm(l + z)^ + Ha) , and if both galaxies have 


-1 
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Vc < 600km/sec, but equals zero otherwise. 

The quantity is plotted in Figure 3 in the case where = lO^K at all times 

(dashed lines) and in the cases I-IV (solid lines). The hrst case closely resembles the one 
considered by Menou et al. (2001) in which they had BHs in all halos (althongh they did 
not apply a detection limit, and hence obtained much higher event rates). The integrated 
event rate IVmrg is listed in each panel of Fignre 3. We hnd that iV^rg is sensitive to the 
reionization redshift, with values of iV^rg ~ 350yr“^ and iV^rg ~ 145yr“^ assnming 2:reion = 7 
and Zj-eion = 12 respectively (cases I and II). As is evident from the plot, the sonrce connts 
may probe the reionization epoch out to large redshifts. There is a drop in the event rate at or 
soon after the reionization redshift. This drop arises because even though some halos below 
the critical temperature still contain BHs until 2 : ~ 2 (2:reion = 7) and 2 : ~ 5 (2:reion = 12) 
respectively, the event rate scales as the sqnare of the BH occupation fraction. Similar 
behavior is seen in case IV, however the smaller valne of ag = 0.6 suppresses the growth 
of structure at early times relative to our hducial case of ag = 0.87. The resulting connts 
are much lower than in cases I-III, N^rg ~ 15yr“^. The drop following reionization is 
much larger if the suppression of gas infall extends to larger galaxies. For example, if 
TTnin = 2.5 X lO^K following reionization, the suppression of the number counts is more than 
an order of magnitude within 2 redshift nnits. 

Similar results are found with Fin, = Face = TArg = 0.01 (JV^rg ~ 340yr“^ and 
A^mrg ~ 140yr“^ for cases I and II). Since mergers quickly result in massive BHs be¬ 
ing nbiqnitous in larger galaxies even if the ocenpation fraction is initially small (Menou, 
Haiman & Narayanan 2001), the nnmber counts and the corresponding spectra (see snb- 
sequent sections) are weakly dependent on Fjni, Face and Fy^rg- We therefore only consider 
Fini = Face = -^mrg = 1 for cases I, H and IV in the remainder of this paper. 

5. The Distribution of Event Durations of Gravitational Waves from 

Inspiraling BH Binaries 

In the previous section we calculated the event rate of detectable gravitational radiation 
sonrees for the LISA satellite. However we would also like to know how long these events 
will last. To this end we calculate the distribution of event durations At for sources having 
observed strains between he and he + Ahc at fc = 10“^Hz. We dehne the event dnration 
At as the time taken for the freqnency to pass throngh a single decade in frequency from 

= IQ-^-^Hz to fe = 10-2-5Hz. 

To calculate the event dnration distribution, we assume that a fraction e^rg of binary 
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BHs produced through galaxy mergers evolve rapidly into the regime where energy loss is 
dominated by gravitational radiation. We can use the merger rate and the mass function of 
halos containing massive BHs together with the BH - halo mass relation in equation (4) to 
hnd the observed number of events of duration At per unit strain per year. Thus 


dhcdAtdt Jq Jq Jo dM dAMdt j^^\dnps/dAMj 1 + z dzdil 

xS [ir(Mbh, AMbh, ;^)] 5 [T(Mbh, AMbh, ;^)] x S{z, Mbh, AMbh), (12) 


where S{x) is the Dirac delta function, from equation (8) 


A(AMbh) = he — S 


/2y/^ G5/3^bhAMbh . .2/3 

VlSy' c^R{z){Moi. + AMbh) 1/3 ^ ’ 


(13) 


and the function S, which was introduced in the previous section controls the fraction of 
galaxy mergers that can result in BH coalescence within a Hubble time. From equation (6) 
we hnd that there is a single redshift corresponding to each delay for a given combination of 
merging BH masses, hence 


T{z) = At- 


40 

’2567r2 


' 2 ' 
,15 


0.5 


C 

hr 


(10-^ (10-^-^ Hz) 


r,-2.5 


Making the assumption that the timescale for the hnal inspiral is short compared with 
a Hubble time (in the emrg x 5* of cases where the satellite sinks to the center of the merger 
product and the resulting BH binary enters the GW regime in less than a Hubble time), 
we can take the mass-function and merger rate terms outside the redshift integral in equa¬ 
tion (12). If we then integrate over 6 [K{M, AM, z)] using J dAM, and over 6 [T{M, AM, z)] 
using J dz we hnd 


d^N 

dhcdAtdt 



d'^hh d N uirg 

dM dAMdt 


[AMbh = K-\0)] 

M 


X 


dK 

dAiffbh 


[AMbh = A-i(O); 


dAMbh \ 

dAM ) 


/ dnhh/dAM\ 

\ drips/dAM ) 


emrg d V 

X dvr-2 ——— 

1 + z dzdil 


X S'(2:, Mbh, AMbh) 



dz 


'z = T-i(O)' 


0[Mbh-A-i(O)], 


(15) 


where we have used the relation 

f dx5[g{x)] = f dx 


^[x = g-\0)] 


-6 bl 


I dx 


( 16 ) 
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The terms are evaluated using the value of AMbh = A“^(0) found by solving the cubic equa¬ 
tion (8), which has only one positive real root. The value = T“^(0) is found numerically 
from equation (14). The Heaviside step function 0 [Mbh — A“^(0)] ensures that mergers are 
not counted twice. Equation (15) has the following intuitive explanation. For each combi¬ 
nation of M and AM there is a merger rate to which the event rate is directly related. The 

derivative ^ dAMbh ) relates the interval of strain to the interval of AM, thus providing 

the number density of mergers ^^^^dAM producing events in dhc- Finally the derivative 
(^) ^ relates the interval of event duration to the interval of redshift in which the source 
can be located. The volume in which these sources will be found is then Air-^^dz. 

dzdil 

Curves representing the number of sources per log At per log he per year are plotted 
in Figure 4 as functions of At for different values of he in cases I-IV. Also shown is the 
confusion limit for observations with no directional information (3.15 x lO^sec/At). Note 
that this confusion line assumes that the detector has no angular resolution; in the case of 
events with locations measured to a positional accuracy A6, the confusion line should be 
lifted by a factor (47r)/(7rAd^). Events with values of strain corresponding to curves below 
this line will be distinct, while those above the line will overlap. Events with lower frequencies 
come preferentially from higher redshift. These longer events are more common, which is 
consistent with the Endings in Figure 3. Events with values of strain he 10-20 - 10-21 
and above will be separable. The figure suggests that in cases I-III around 10 events per 
year with he ~ 10-2° will be observed, with each lasting around a few days. At a larger 
amplitude, he ~ lO-i®, about one event per year will be observed lasting less than a day. 
Even larger amplitude events from inspirals will be very rare. In case IV we find that the 
lower event rate means that only a few distinct events are expected per year at strains of 
he ~ 10-20 of more. Another prominent feature of case IV is the broadness of the curve for 
he = 10-20, being dominated by sources producing the peak in number counts near z ^ 4 
(see Figure 3). 


6. The Spectrum of the Gravitational Wave Background from Inspiraling BH 

Binaries 

LISA will be sensitive to the final in-spirals and/or ringdown phases of most massive BH 
mergers (Haehnelt 1994; Hughes et al. 2001). However, with the exception of the occupation 
fraction, the event rate does not directly carry information on properties of the BH population 
itself such as the slope and normalization of the BH mass-function. Individual events will 
yield useful information (Hughes 2002), but we might also expect the number counts of 
sources in frequency and strain (i.e. the time averaged number of sources per unit frequency 
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per unit strain at a single epoch) to provide constraints on the relation between BH and halo 
mass. 


In analogy with the previous section we can use the merger rate and the mass function 
of halos containing massive BHs together with the BH - halo mass relation to hnd the 
instantaneous, time averaged number counts that should be observed in a snapshot of the 
sky per unit frequency per unit strain per co-moving Mpc, 


dhcdfc 


z = 


dM / dAM 


, dnhh d^N„ 


dM dAMdt 


X S'(/c, 2:,Mbh, AMbh) 


/ dn^h/dAM 

\dnps/dAM 
where from equation (6) 


5 [ir(Mbh, AMbh, t - Wt)] 5 [i^(Mbh, AMbh, t - 4tart)] , (17) 


F{t) ^ /e - 


(1 + z) 2 tt 


,_8/3 256 MbhAMbh ,, , J 

^ 5 c5 Mbh + AMbh 


Three conditions are imposed on the number counts and are incorporated in the function 
S. The hrst requirement is that the circular velocity of the merger product be smaller than 
600km/sec. In § 4 we assumed that coalescence and the emission of GWs at fc = 10“^Hz 
would only occur where the merging secondary galaxy sinks before the next comparable 
merger event and in less than a Hubble time. For calculation of equation (17) we note that 
within the GW dominated regime there are some frequencies and mass combinations for 
which coalescence does not occur within a Hubble time. We therefore assume that a BH 
binary produced through a galaxy merger will only result in the emission of gravitational 
radiation at a frequency fc if the binary coalescence as well as the prior sinking of the satellite 
occur before the next comparable merger event and within a Hubble time. Thus the function 

S{z, Mbh, AMbh) is set to max (^AM [fdecay + f(/c) - f(/max)] ciAMdt ^ satellite 

galaxy’s orbital decay time tdecay plus the inspiral time [t(/c) — f(/max)] is smaller than the 
Hubble time, and if both galaxies have Vc < 600km/sec, but equals zero otherwise. 


As before we make the assumption that the timescale for the hnal inspiral is short 
compared with a Hubble time and take the mass-function and merger rate terms outside 
the time integral in equation (17). If we then integrate over 6 [K{M, AM, f — Gtart)] using 
f dAM, and then over 6 [^(M, AM, t — Gtart)] using f dt we hnd 




log^o[d^N/dlog(h^)dlog(f^)dz] log^Q[d^N/dlog(h^)dlog(f^)dz] 




logio(^) 


logio(hc) 


Zreion=12 T^in = 250K F=0.0 ag = 0.87 


Zreion=12 T^i„=10000K F=1.0 ag= 



logio(hc) 



Fig. 5. — Plots of diog{h)diog{f ) different frequencies. In each panel the thin lines (from top to b' 
correspond to fc = 10“^°, 10“®, 10“®, 10“^, 10“®, while the thick lines correspond to 10“^, lO""* 
10“^ and 10“^Hz. The panels show cases TIV whose primary parameters are listed at the top of eac 
We assume cq = 10“^ '^ and 7 = 5. 

































- 21 - 


The derivative ^[t = F ^(0)] is easily computed by substituting fc and AMbh = K ^(0) 
into which is in turn obtained from equation (6). The Heaviside step function 

0 [Mbh — -^'“^(O)] ensures that mergers are not counted twice. 

As before, equation (19) has the following intuitive explanation. For each combination 
of M and AM the derivative ^ dAMhh ) rslates the interval of strain to the interval of AM, 

thus providing the number of mergers per volume per time, ^^^^dAM, producing events 
within dhc- The derivative (^) ^ relates the interval of frequency to the time during which 
the binary is in that frequency interval. The number count in a snapshot of the sky is 
proportional to the number of mergers during that time. 

The number counts of GW sources per logarithm of frequency per logarithm of strain 
over the sky is 

dl0g(he)dl0g(/c) “ Jo ^ ^ ^ 

This quantity is plotted in Figure 5 as a function of he for different values of fc and for cases 
I-IV. As expected, the number counts increase toward low frequencies since the inspiral 
lifetime scales as t oc However at all frequencies the number counts are cut off at 

both high and low amplitudes. The reasons are as follows. First, the cutoff at large he is 
caused by the fact that the largest amplitude waves are produced at the hnal phase of the 
inspirals. Since the inspiral has a maximum frequency (Hughes 2002), reached when the 
binary separation becomes comparable to the Schwarzschild radius, each binary also has a 
corresponding maximum amplitude. This effect plays a significant role in the characteristic 
strain spectrum in the frequency range applicable to LISA, a point to which we will return 
in the following section, § 6.1. Second, the cutoff at low amplitude is caused by the fact 
that at a fixed frequency, he oc MbhAMbh(Mbh + AMbh)”^'^^. Since both Mbh and AMbh 
have a minimum value Mbh,min, the strain must be larger than a value proportional to 


6.1. The Characteristic Strain Spectrum 

The characteristic strain spectrum, hspec[/c] (e.g. Phinney 2001; Jaffe & Backer 2002), 
describes the spatial spectrum of the stochastic background of gravitational radiation. The 
characteristic strain spectrum is related to the strain power-spectrum Ah through 


f'spec(fc) = VfeShife)- 


( 21 ) 
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where 5'h(/c) is calculated from 


Stifc) 




dz hi 


' dhcdfc 


{z)4:TI 


d^V 

dfldz 


( 22 ) 


The contribution to the characteristic strain spectrum from sources in a redshift interval can 
be calculated by integrating between the desired limits. 

We have plotted hspedfc) in Figure 6 for cases I-IV. The thin lines from bottom to top 
show hspecifc) in the redshift intervals z > 6, 2 < z < 6, and z < 2. The solid thick lines 
show the total hspec(/c)- The characteristic strain spectrum is dominated by sources having 
z < 2. The thick dashed line shows the spectrum deduced by Jaffe & Backer (2002). They 
found hspecifc) 10 for the region accessible to pulsar timing experiments. Figure 6 

shows their result at nHz frequencies and extrapolates it to the higher frequencies probed 
by LISA for comparison. 

At frequencies between 10“®Hz and 10“®Hz, the characteristic strain spectrum has a 
logarithmic slope of -2/3. This slope results from the dependence of energy loss on fre¬ 
quency as explained by Phinney (2001). At the higher frequencies relevant for LISA, we 
hnd departures from the -2/3 slope. The steepening of hspedfc) is due to the loss of power 
at high frequencies associated with the maximum frequency during inspiral. At frequencies 
below 10“®Hz there is loss of power arising because many very low frequency sources do not 
coalesce within a Hubble time. The thick grey lines in the lower right corner of each panel 
in Figure 6 show the instrumental noise equivalent of hgpec computed from the instrumental 
noise spectral amplitude S'“®*(/c). The curves are from a fit (Hughes 2002) to calculations 
of *S'™®*(/c) (Folkner 1998). The amplitude of the characteristic strain spectrum should be 
detectable by LISA over the frequency range ~ 10“^-10“^Hz. However the low numbers 
of events will make the power-spectrum noisy in this regime. In general the binaries can 
be resolved provided less than 20% of the frequency bins are hlled at any one time. As we 
showed in § 5, most events are of sufficiently short duration to allow each to be resolved sep¬ 
arately. The white-dwarf/white-dwarf binary foreground was estimated by Phinney (2001) 
and is plotted as the thick dotted lines [cut off below 10“‘^Hz since low frequency binaries are 
not expected to contribute significant gravitational waves (Hils, Bender, & Webbink 1991). 
Through most of the LISA band, this foreground is more than an order of magnitude below 
the expected background due to coalescing massive BH binaries. 

We also show in Figure 6 the current upper limits on the GW background (large dots). 
The most recent limits are quoted by Lommen (2002). They find a limit on the fractional 
energy density in gravitational waves per logarithm of frequency at /c = 1.9 x 10“®Hz is 
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Zreion = 7 T^i,= 10000K F=1.0 as = 0.87 z^eion=12 T^i,= 10000K F=1.0 ag = 0.87 



fc (Hz) (Hz) 


Fig. 6.— The characteristic strain spectrum. In each panel the thin lines in order from bottom to top 
show the spectrum in different redshift intervals and the thick dark solid line is the total strain spectrum. 
The thick dashed line shows the characteristic spectrum hspec = found by Jaffe & Backer (2002), 

and the thick dotted line the estimated stellar binary foreground from Phinney (2001). The thick grey solid 
line in the lower right corner of each panel shows the expected instrumental noise for LISA. The dot shows 
the current limit from pulsar timing. The panels show cases I-IV whose primary parameters are listed at 
the top of each plot. We assume eo = 10“^'^ and 7 = 5. 
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f2g/i^ < 2 X 10“®. We then find a corresponding limit^ on the characteristic strain spectrum 
of hspec(1.9 X 10“®Hz) < 3 X 10“^^. The limits are close to the theoretical spectrum in cases 
I-III. We will return to this important and unexpected result below. 

Our results in cases I-III for the amplitude of the characteristic strain spectrum are 
different from the findings of Jaffe & Backer (2002) in the regime accessible to pulsar timing. 
In particular we find that the models predict hspec(l-9 x 10“®Hz) that is an order of magnitude 
higher than predicted by Jaffe & Backer (2002), and which is close to the limits of detection 
by the pulsar timing measurements. There are several possible causes for this difference. Jaffe 
& Backer (2002) use a power-law extrapolation of an estimate of the spheroid merger rate 
to higher redshifts, while in fact galaxies were made of smaller sub-units that had a higher 
merger probability at earlier cosmic times. Also, they use a phenomenological prescription 
for the merger rate that may substantially underestimate the actual merger rate for short 
duty cycle events. Furthermore, they assume the same merger rate for all galaxies, and 
do not assign a relation between black-hole and galaxy mass. In essence, we have assumed 
BHs to be associated with dark-matter halos, for which the merger history is understood 
rather than use an estimate of the average spheroid merger rate. The two assumptions lead 
to different results because the total merger rate of galaxies above the cooling threshhold 
as computed from the Press-Schechter formalism is higher than estimates of the observed 
spheroid merger rate. We find that if we artificially force equality between the merger rates 
computed in our work and those used in Jaffe & Backer (2002), we obtain a similar spectrum 
well below current limits. 

The results should also be dependent on the BH - halo mass relation. The amplitude 
of hgpec is proportional to [see equation ( 8 ), ( 21 ) and ( 22 )]. As a result, a lower value 
of €o will result in hspec having a lower amplitude. Moreover, the production of power at 
very low frequencies is dominated by very large BHs, and so hspec should be sensitive to 7 , 
with larger 7 values resulting in larger amplitudes. We note that the log-normal distribution 
employed by Jaffe & Backer (2002) is steeper at the highest masses than our BH mass 
function computed from the Press-Schechter mass function and the Mbh — M relation, which 
might further contribute to our differing results. We will return to the dependence on the 
BH - halo mass relation in the following section (§ 6.2). 

Our simulation of Case IV yields a value of hspec(l-9 x 10“®Hz) that is an order of 
magnitude below current limits, and which is more consistent with the findings of Jaffe & 
Backer (2002). Thus, the power in GWs is quite sensitive to the value of (Jq. Indeed hgpec is 


®This value is an order of magnitude larger than that plotted in Figure 8 of Jaffe & Backer (2002). Use 
of the value in their paper leads to even tighter limits on the parameters ep, 7 and Cmrg 
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significantly more sensitive to ag than to the details of the merger history provided that the 
history results in BHs being ubiquitous in massive galaxies in the local universe. 

We hud that the characteristic strain spectrum at low frequencies is dominated by 
galaxies having masses larger than ~ merging at redshifts below the peak in quasar 

activity {z < 2). This is particularly true of case IV. Theoretical studies of the quasar 
luminosity function hnd that the lack of cold gas at low redshifts explains the drop in quasar 
activity in massive galaxies® (Kauffmann & Haehnelt 2000). If as suggested (e.g. Gould & 
Rix 2000) gas that is driven to the center of the galaxy during a merger is the mechanism 
that brings the BH binary into the GW dominated regime, then we might expect galaxy 
mergers to be inefficient (e^rg 1) at producing gravitational radiation at low redshifts. If 
emrg "C 1 , then we expect the amplitude of hgpec which is proportional to ^emrg, to be lower 
than predicted. Our results from cases I-III imply that limits on the efficiency can already 
be placed for some values of 7 and eo using limits on hspec from pulsar timing experiments 
(§ 6.2). If emrg is found to have a value of order unity through the detection of the background 
by future pulsar timing experiments or by the LISA satellite, then this would strengthen the 
case for similar values at yet higher redshifts. The number counts of mergers detected by 
LISA will therefore directly determine the BH occupation fraction at high redshift. 


6.2. Limits on eo and 7 

The models considered in this paper predict characteristic strain spectra with values 
at /c = 1.9 X 10“®Hz that are close to the current limits from pulsar timing experiments 
(Kaspi et ah 1994; Thorsett & Dewey 1996 and Lommen 2002). So far we have concentrated 
on the spectrum and event rate for eo = 10 “®'^ and 7 = 5 , which are the values favored 
by simple feedback arguments and the theoretical model for the quasar luminosity function 
described in Wyithe & Loeb (2002). However variation of eo and 7 has a signihcant effect 
on the GW signal. Since hspec is dominated by the mergers of massive (> IO^^Mq) galaxies 
at low redshifts {z < 2) in which supermassive BHs are thought to be ubiquitous, it should 
be quite insensitive to the details of the uncertain evolution in the BH mass function at 
higher redshifts. The amplitude of hspec(l-9 x 10“®Hz) therefore provides robust constraints 
on combinations of the parameters emrg, 7 and eo. 


®It is worth noting that the fiducial values of eo and 7 used in this work are taken from results of the 
theoretical quasar luminosity function described in Wyithe & Loeb (2002) which over-predicts the number 
of bright quasars at low redshift. However this approach does not over-predict the GW back-ground by 
predicting the existence of to many very massive BHs because similar values of eo and 7 have been measured 
by Ferrarese (2000) in nearby galaxies. 
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Fig. 7.— Contours of the value of the characteristic strain spectrum at /c = 1.9 x 10 “®Hz as a function 
of €o and 7 . The thick line shows the current upper limit from pulsar timing experiments. The dot shows 
the position of the parameter pair eg = 10“®'^, 7 = 6 . The open circle the values found by Ferrarese (2002) 
assuming that the circular velocity equals the virial velocity (7 = 5.46, eg = 10“^ °®), while the square shows 
the observed values (7 = 4.95, cq = 10“'^'®"') assuming an NFW profile (Bullock et al. 2001). The panels 
show cases TIV whose primary parameters are listed at the top of each plot. 
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In Figure 7 we show contours of hspec(l-9 x 10“®Hz) as a function of eo and 7. The thick 
grey line shows the contour at the current limit of hgpec < 3 x 10“^^. The dot shows the 
values Cq = 10“^'^ and 7 = 5. As mentioned the value is close to the observational limits. 
The open circle shows the values found by Ferrarese (2002) assuming the circular velocity 
equals the virial velocity (7 = 5.46, eo = 10“®'°®), which is the assumption implicit within 
our formalism. The square shows the corresponding observed values (7 = 4.95, eo = 10“^'^"^) 
obtained assuming a Navarro, Frenk & White (1997) prohle (Bullock et ah 2001). We see 
that if coalescence is prompt for all binaries, then the pulsar measurements already exclude 
parameter space close to the currently favorable values. 

Adoption of the values determined by Ferrarese (2002) for the BH - halo mass relation 
implies either that the background will soon be detected, or that the efficiency for BH binary 
coalescence at low redshift is emrg < 1- Yu (2002), investigated the evolution of hypothetical 
BH binaries in realistic models of nearby early-type galaxies. She found that BH binaries 
are likely to survive for longer than a Hubble time in spherical or weakly flattened/triaxial 
galaxies. No evidence has yet been found for BH binaries in the centers of the galaxies 
studied. The upper limit of the range of semimajor axes of surviving binary black holes 
as determined by Yu (2002) is close to the HST resolution for typical nearby galaxies (i.e., 
galaxies at Virgo). Thus, the absence of evidence for double nuclei in nearby galactic centers 
does not necessarily mean that they have no BH binaries (Q. Yu, private communication). 
Further improvement by an order of magnitude in the upper limits on the GW background 
will imply that emrg < 1 and that BH binaries may survive in galaxies rather than coalesce. 
Yu (2002) also computed the periods of surviving BH binaries. Of particular importance to 
this study is that nearly all of the surviving BH binaries would have periods longer than 20 
yrs, and would therefore not contribute to hspec(l-9 x Hz). 


7. Conclusions 

We have computed the characteristic strain spectrum for GWs, as well as the number of 
events detectable by the planned LISA observatory due to coalescing massive BH binaries. 
Rather than adopt a phenomenological approach which is limited by the redshift horizon of 
current observations, we calculated the expected merger history of galaxies out to arbitrarily 
high redshifts under different assumptions. Our model uses the observed relation (Ferrarese 
2002) between the circular velocity of galaxies and their central BH mass, and assumes that 
massive BHs form in galaxies into which gas can condense and cool. An important ingredient 
in the prescription is the observation that the slope and normalization of the relation holds 
out to high redshift (Shields et al. 2002). The underlying principles of this approach have 



been shown to reproduce current data on the luminosity function of quasars (Wyithe & Loeb 
2002; Volonteri, Haardt & Madau 2002). 

We hud that if BH formation is ongoing in galaxies with masses above the cooling 
threshhold, then reionization produces a sharp drop in the number counts of LISA sources 
as a function of source redshift (see Fig. 3). In the instance that all binary BHs proceed to 
coalescence our models generically predict that hundreds of events with strains he > 10 “^^ 
will be detectable per year by LISA. The strongest sources have short durations and will not 
suffer from a confusion limit. Interestingly, most of these sources originate from redshifts 
z>7. 

In the uHz regime accessible to pulsar timing experiments, we hnd that the amplitude 
of the characteristic strain spectrum (related to the power in the gravitational wave back¬ 
ground) is close to current limits assuming values (favored by Wyithe & Loeb 2002) of 7 = 5 
and eo = 10“®'^ for the slope and normalization relating the BH to host galactic halo mass 
(see equation 4). Based on local observational data on early type galaxies, Ferrarese (2002) 
hnds 7 = 5.46 and Cq = 10“^-°® (where the circular velocity is assumed to equal the virial 
velocity), or 7 = 4.95 and Cq = 10“^'^"^ (where the halo density is assumed to follow an NFW 
prohle). Calculations of the characteristic strain spectrum assuming these parameter sets 
are comparable to current limits, and already rule out only slightly different combinations of 
7 and eo if all binary black holes proceed to coalescence. Alternatively, future improvements 
in the observational limits will allow limits to be placed on the fraction of black hole bina¬ 
ries that proceed to coalescence from calculations of the strain spectrum that use measured 
values of 7 and cq. 

We hnd that the amplitude of the characteristic spectrum of strain is dominated by the 
mergers of massive galaxies (M > IO^^Mq) at low redshift {z < 2). As a result, hgpec is 
quite sensitive to the value of Since BHs are ubiquitous in the local universe, detection 
of the background either by future pulsar timing measurements or by the LISA satellite will 
determine the efficiency with which binary BHs coalesce. 

Observations of GWs are likely to provide insight into the formation history of massive 
BHs. LISA will open a new window to studies of the high redshift universe. The number 
counts of gravitational wave sources will provide an unprecedented probe of the reionization 
epoch, and constrain the merger rate of galaxies and the growth history of black holes at 
early epochs that are not yet probed by conventional astronomical telescopes. Pulsar timing 
arrays are suitable to measure the background arising from BH binaries at lower redshifts. 
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